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ABSTRACT 

A two-way time synchronization experiment was per- 
formed in the spring of 1979 and I 98 O via the Ital- 
ian SIRIO-1 experimental telecomiminications satel- 
lite. 

The experiment was designed and implemented by the 
Istituto Elettrotecnico Nazionale, Torino (Italy), 
to precisely monitor the satellite motion and to 
evaluate the possibility of performing a high-preci- 
sion, two-way time synchronization using a single 
communication channel, time-i^ared between the par- 
ticipating sites. 

The results of the experiment show that the preci- 
sion of the time synchronization is between 1 and 
5 ns, while the evaluation and correction of the 
satellite motion effect has been performed with an 
accuracy of a few nanoseconds or better over a time 
interval from 1 up to 20 seconds 


INTRODUCTION 

The principal features of the SIRIO-1 time synchronization 
experiment can be briefly summarized as follows; 

- the experiment was designed to precisely monitor the satel- 
lite motion and the effects of this motion on the time syn- 
chronization accuracy; 
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- the time synchronization is performed using the two-way 
time synchronization technique and a single communication 
channel, time-shared between the two sites; 

- the experiment tests a new technique that has been proposed 
to correct for the satellite motion effect while performing 
the time synchronization; 

- by using a single communication channel, no effects affect- 
ing the accuracy of the tim^ synchronization at the 1 ns 
level are due to the space segment (from one ground anten- 
na to the other) , thanks to the high frequencies used for 
the RF carriers; 

- the time signals used allow the independent determination 
of the uncertainties of the time-of-arrival measurements 
at the two stations, to separate the contribution of each 
station to the total precision. 

This last feature can be important to understand the contri- 
bution of local phenomena (ground equipment, atmospheric con- 
ditions affecting the signal attenuation, especially rain, 
etc.) to the synchronization precision. 

ORGANIZATION OF THE EXPERIMENT 

A detailed description is given in ref. 1. Only a few remarks 
are given here, mainly for reference purposes. 

Two groiHid stations, Fucino and Lario (fig. 1), participate 
to the experiment. Both sites are in Italy, in the north- 
ern (Lario) and in the central part (Fucino) of the country; 
one lEN Cesiiira clock was installed at each site. 

Fucino (fig. 2) transmits its time signal at 0 seconds of the 
synchronization frame (Fucino time), acting as station ’A', 
while Lario transmits its own signal at 0.5 s (Lario time), 
acting as station *B* , 

Two times-of-reception are then measured at each site; with 
reference to fig. 3, these are T^ and T^ at Fucino and T^ 
and T^ at Lario (we have no need to measure Tq and T^ since 
these are, known). 
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To simplify the notation and for a better understanding, we 
may note that T4 - Tq is the time of propagation of the time 
signal from Fucino to Fucino, T-| - Tq is (neglecting for now 
e) the time of propagation from Fucino to Lario, etc., so we 
can write; 

T(FF) = T4 - T-, = T4 
T(FL) T^ -Tq^^ Tt 

^ ^ T(LF) S' T3 - T2**»i T.- 0.5 

T(LL) = T5-T2 = T5-O.5 

(*) local time reference 

T(FF), T(FL), T(LF) and T(LL) are the actual results of the 
time measurements at the two sites [ except for the subtrac- 
tion of 0.1 s, resulting from the hardware implementation, 
see ref. 1 ] and will be used in any following computation. 

Two data types are considered: pseudo-range data, such as 
T(FF) and T(LL), that are the time intervals measured against 
the same time reference, and synchronization data, such as 
T(LF) and T(FL); the starting and ending times of the latter 
intervals are measured with reference to different clocks. 

Data format 

Actually, each one of the values listed in ( 1 ) results from 
the measurement process as the mean over ten independent 
measurements: a rough data file is shown in fig. 5 . Each 
time of arrival is then evaluated as the arithmetic mean of 
the measured data. The data is rejected if the associated 
standard deviation is larger than 100 ns; however, less than 
0,5f° of the data was rejected because their standard devia- 
tion was exceedingly large. On the average, the standard de- 
viation for each of the tirae-of-arrival evaluations, based 
on ten data values, is in the range 10 to 50 ns. 

The basic synchronization frame lasts 1 s and is repeated 
every 10 s; during the 1979 series of measurements, small 
groups of data (10 to 15 measurement frames) were recorded 
sequentially, to characterize the performance of this tech- 
nique over time intervals of 100 to 150 seconds: this was 
actually performed also to verify the assumption of a linear 
motion of the satellite and the validity of the correction 
used (see eq. (7) to (15), ref. 1) over this time interval. 
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During 1980 a second series of measurements were per- 
formed over longer time intervals (up to 16 minutes); the 
comparison between the results obtained via the satellite, 
the TV method and portable clock trips are shown in fig. 6, 
where c ^ is the difference between the two atomic clocks 
located at the groimd sites. An expanded view over four con- 
secutive days of satellite measurements is given in fig. 7. 

The overall accuracy of the clocks comparison was estimated 
to be between 50 and 100 ns, with reference to some portable 
clock comparisons. 


Synchronization estimate (x = 1 s) 

The clocks difference e-fc time t, defined as; 


(2) e. = t(B) - t(A) = t(LAR) - t(PUC) 

X 


is given (see ref. 3), over the basic synchronization frame 
time interval! (x = 1 s), by the equation: 

T (PL) - T (LP) 

(3) e.(x = 1 s) = — i 2 + (0.5)*C 

^ 2 
since (tg-t^ ) is 0.5 s (see fig. 3). 


The range-rate correction C to e^, as defined in ref. 1, 
computed as: 


(4) 


C 



20 


is 


The magnitude of C was usually found to be in the range 
2*4 ns/s, yielding a range rate correction of 1 to 2 ns over 
the basic frame. 


An estimate of e is obtained by taking the arithmetic mean 
of a number of successive data (from 100 s up to 16 minutes). 
If e would be constant over this measurement interval, then 
the standard deviation <r(e ) of the data wo\ild be related to 
the precision of the method. 

The evaluation of c and ^(e ) is carried on by applying a 
3-sigma width filter; that is, the £ -j; value is rejected if 
the residual | e-e . |> 3cr ; if anye . has been rejected, a new 

X X 
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e and O' (e) are evaluated using the remaining data. 

This procedure is repeated \mtil no more data are rejected; 
usually this filter rejects less than the 5/» of the available 
data. 


Assuming a nonnal distribution of the data, the one-sided 
estimated error 1 6c | in the determination of c is given, with 
the 99. confidence, as; 


(5) 


5cl = 


‘^•'^ 0 . 995 '''-^’ 




The magnitude of 6c , when c is computed over 50 to 100 data 
c-t, ranges usually between 1 and 5 ns (see fig. 7 and fig. 8). 

The capability to look at the precision of the method over 
short time intervals is demonstrated by fig. 8 and fig. 9, 
where two-hours data are plotted, together with an \mweighted 
least- squares fit of the available data. 


As it can be seen, the average error of the fit is quite small, 
less than 1 .5 ns. 


Expanded time synchronization frame (x = 10 s and X = 20 s) 

As explained in ref. 1 , the time synchronization frame can 
be lenghtened by simply rearranging the data. 

This shows the capability of the method to synchronize two 
clocks by using a single commimication channel even if a fast 
switching of the RF carrier at the two sites is not possible 

and a large effect due to the satellite motion is then expect- 
ed; however, the amount of the correction due to this motion 

can be computed very accurately by using the pseudo-range 
data available. 


In this case, the clocks difference can be computed as: 

T (FL) - T (LF) 

(6) e^(x = 10 s) = + (10.5)*C 


(since now (t 2 -t-| ) a: 10 .5 s) where C is given by eq . (4). 

This is equivalent to have a station transmitting its time 
signal at to and the other site transmitting at (to+10.5)s. 
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Accordingly, to simulate a longer time interval (x = 20 s) we 
can write: 


e^(x = 20 s) = 




+ ( 20 . 5)-0 


where now C is computed as: 


(8) 


C 


1 

2 40 


A comparison of single measurements of e^, evaluated by using 
eq . (3)f (6) and (7) is presented in fig. 10, over a 200 s 
time interval. The agreement between the values of for dif- 
ferent x's is remarkable, if we note that the correction to 
be applied to eq. (6) and (7) amoixnts respectively to about 
48 and 95 ns in most cases. 

Further analysis on the experimental data 

Instead of using the differences between the measurement data 
T(FL) and T(LF) to compute e-j;, it is possible to use a poly- 
nomial fit of the data and then evaluate over the coef- 
ficients of the fitted polynomial. 

This procedure has the advantage that less data are exchanged 
between the two sites (the fit coefficients only) and that 
also missing data points at one site can be recovered from 
the fitted curve. 

By writing: 

T(FF) = f. a. (t-to)^ 

T(LL) = f. b. (t-to)^ 

(9) g . 

T(FL) = Zi c.(t-to)^ 

T(LF) = I. d (t-to)^ 

q-L 1 

eq . (3) can be written as: 

(10) e^(x = 1 s) = ^ ^Zi(c^-d^)(t-tQ)^+(0.5) *C 

where: 
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( 11 ) 


c 



^ i ( a^+b^ ) ( t-tQ^ 1 0 ) “ oi ^ ^ ^ ^ ^ 

20 

fi (a.-Hb.)«10^ 

20 


If is constant over the measurement interval, 
and eq. (10) becomes: 


then e 


t ^to 


(12) 1 (dQ-c^) + (0.5)*C 

If only the linear (velocity) terms are significant and any 
higher order term (acceleration) of the satellite motion is 
neglected, then 0 is given by 

a. + b. 

( 13 ) 0 = + ^ ^ 


This procedure was actually carried on over measurements in- 
tervals up to 600 s; over time intervals up to 200 s it was 
fovind that a linear (first order) fit (as given by eq. (12) 
and ( 13 )) was usually good anough to evaluate z. and any 
further increase in the degree of the polynomial does not 
improve the fit; obviously this result depends mainly on the 
clocks and the synchronization process behaviour. 

This was also a further check of the correctness of the lin- 
ear motion assumption as given by eq . (4). 

The test of statistical significance (see ref. 2) of the com- 
puted coefficients of the polynomials (9) was carried on by 
computing the standard deviation of the estimated coeffi- 
cients. 


This was done in the following way: the least squares fit is 
expressed by the normal equations that, in matrix form, are: 


(14) (X'X)(>> = (X'Y) 

where: (X'X) is the normal equations matrix (or x-products 

matrix); 
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(X'Y) is the cross-products matrix; 

(p is the vector of the coefficients. 


An estimated value of the standard deviation of the fit is 
computed as: 


(15) 


2 

O' 


1 

(n-k) 


2 


r.r. 

1 1 


where n is the niimber of data, k is the number of estimated 
coefficients and r^^ is the i-th residual. The estimated stan- 
dard deviation (ref. 2) of the k-th coefficient is given by: 

“’k = ’’■T^ 

_1 

where c^j^ is the k-th diagonal element of the matrix (X*X) . 
Actually the computation of the (X*X)“^ matrix was performed 
with Gauss- Jordan reduction and pivot search to minimize the 
numerical computation errors. Then an estimated confidence 
interval for the coefficient can be computed, by using the 
Student t-distribution at (n-k) degrees of freedom. 


Ground-equipment delays measurements 

Two types of measurements were performed: test-loop-translator 
(TLT) measurements and transmitting-chain delay measurements 
(LARIO site only) . 


Test-loop-translator measurements 

The experimental set-up is shown in fig. 11. The measurements 
performed showed a precision around 1 ns and a long term (1 
month) stability of the delay in the order of 1 to 3 ns, if 
the groxmd equipment is operated at the same power level. The 
measurements were performed in the same operating conditions 
as during the synchronization sessions. 

The loop-delay was foiuid to be 3.776 ;us (LARIO) and 4.618 ;is 
(PUCINO) on the average. 

Transmitter delay measurements (LARIO site only) 

The proposed use of a microwave cavity as a frequency discrim- 
inator was tested. The cavity was characterized by a Q of 
150.0 at the transmission frequency. Unfortunately, the only 
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way to couple the cavity near the TLT was via an existing 
directional coupler, that attenuates the signal more than 
40 dB. 

The rectified signal consequently was very small, since the 
cavity contributed an additional attenuation of about 12 dB, 
and, under these conditions, the measurement was impossible. 

In order to increase the rectified signal, it was necessary 
to increase the frequency deviation: in these conditions the 
communication equipment was working outside the range of nor- 
mal operation and the overall system response degraded no- 
ticeably. The measurement jitter, for instance, increased up 
to 50 ns (1 signa), as compared to the 1 ns found in the TLT 
measurements; this was verified by performing the same TLT 
measurement, but with the larger frequency deviation. 

The reproducibility of the measurements, mainly related to 
the critical setting of the microwave cavity (the cavity res- 
onance was adjusted at the RF carrier frequency when no mod- 
ulation was applied), was better than 100 ns, even when 
working in these very critical conditions. 
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Fig. 3 - Timing diagram 
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Fig. 4 - The LARIO station antenna 
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26.020 


Fig. 7 

SIRIO-1 TIME SYNCHRONIZATION 
DATA 
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25.790 
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e (ps) 
25.892 
.890 
.888 
.886 
.884 
.882 
.880 
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£(> 15 ) 

25.800 

.796 

.792 

.788 

.784 

.780 

.776 

.772 
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2- transmitter delay only 
Fig. 11 - Ground equipment delay measurement set-up 








Fig. 13 - Automatic TV measurements subsystem 
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QUESTIONS AND ANSWERS 


MR. DAVID ALLEN, National Bureau of Standards 

Two questions; one what Is the size of the antenna involved? 
PROFESSOR LESCHICUTTA: 

Yes, please. In the ground station the size is 17 meters with a 
true bandwidth -- RF bandwidth of 34 megahertz. Also, in the re- 
peater satellite in the real bandwidth of base band bandwidth of 
6 megahertz. 

As regards the experiment on the ship, the diameter of the 
dish is on the order of 2 1/2 meters but the bandwidth is just 
1.5 megahertz. So, obviosuly the precision should be deteriorated. 

MR. ALLEN: 

I thought it was a very excellent result that you received. I have 
one question in regard to the equation. Because the stations are 
basically north/ south you would not see any effect due to the SANYAC 
correction. 

PROFESSOR LESCHIUTTA: 

Yes, the SANYAC correction is 13.5 nanoseconds, in our case because 
the area of the path is very small in the equitorial plan. 

MR. ALLEN: 

The correction was not in the equation? 

PROFESSOR LESCHIUHA: 

No, no, it was not included but is in the order of 14 nanoseconds. 
MR. ALLEN: 

Thank you. 

PROFESSOR LESCHIUTTA: 

13 foot. 

MR. ALLEN: 

Very good . 
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CHAIRMAN BUISSON: 

Any other questions? 

MR. LAUREN RUEGER, The Johns Hopkins University/ Applied Physics Laboratory 

Do you ever take advantage of planning your experiments on the satel- 
lite motion when the relative changes to the stations are minimized? 

PROFESSOR LESCHIUTTA: 

There, again, we are not in a position to do so. We just receive 
for some hours during the day, but we are planning periods of ex- 
periments to make all day measurements in order to follow the satel- 
lite. 

In previous experiments we have seen the maximum relative 
speed of the satellite is of the order of 3.5, 4 meter per second. 
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